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ABSTRACT 



We present the results of our broadband spectral analysis of 42 
SGR J1550-5418 bursts simultaneously detected with the Swift /X-ray Telescope 
(XRT) and the Fermi/Gamma-ray Burst Monitor (GBM), during the 2009 January 
active episode of the source. The unique spectral and temporal capabilities of the XRT 
■ Windowed Timing mode have allowed us to extend the GBM spectral coverage for 

these events down to the X-ray domain (0.5 — 10 keV). Our earlier analysis of the GBM 
data found that the SGR J1550-5418 burst spectra were described equally well with a 
Comptonized model or with two blackbody functions; the two models were statistically 
indistinguishable. Our new broadband (0.5 — 200 keV) spectral fits show that, on 
average, the burst spectra are better described with two blackbody functions than with 
the Comptonized model. Thus, our joint XRT/GBM analysis clearly shows for the first 
time that the SGR J1550-5418 burst spectra might naturally be expected to exhibit a 
more truly thermalized character, such as a two-blackbody or even a multi-blackbody 
signal. Using the Swift and RXTE timing ephemeris for SGR J1550-5418 we construct 
the distribution of the XRT burst counts with spin phase and find that it is not 
correlated with the persistent X-ray emission pulse phase from SGR J1550-5418. These 
results indicate that the burst emitting sites on the neutron star need not be co-located 
with hot spots emitting the bulk of the persistent X-ray emission. Finally, we show 
that there is a significant pulse phase dependence of the XRT burst counts, likely 
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demonstrating that the surface magnetic field of SGR J1550-5418 is not uniform over 
the emission zone, since it is anticipated that regions with stronger surface magnetic 
field could trigger bursts more efficiently. 

Subject headings: pulsars: individual (SGR J1550-5418, IE 1547.0 - 5408, PSRJ1550 - 
5418) - stars: neutron - X-rays: bursts 



1. Introduction 



Soft Gamma Repeaters (SGRs) and Anomalous X-ray Pulsars (AXPs) are the observational 
manifestations of magnetars - isolated neutron stars p ossessing extreme magnetic fields, B > 10 G 
(jDuncan k, Thompson! Il992l ; iKouveliotou et al.lll998l ). Besides being bright X-ray sources, SGRs 
and AXPs emit intense bursts in hard X-rays / soft 7-rays on a highly unpredictable frequency with 
peak luminosities ranging from 10 38 erg s _1 to > 10 47 erg s _1 . These energeti c events are attributed 
to th e cracking of the solid neutron star crus t by magnetic s tress build-up ([Thompson Duncan 
1995 ) or t o magnetic field line r econn ection ( Lvutikovj 2003). For detailed reviews on SGRs and 
AXPs, see lWoods fc Thompsonl (l2006h and iMereghettil (boosh . 



SGR J1550-5418 was discovered as a point source with the Ein stein Observatory while search 
ing for X-ray emission from radio emitting supernova remnants (jLamb &; Markert 



1981 



The 

source was suggested as a magnetar candidate by the similarity of its persistent X-ray spec - 
trum to AXPs and its association with a young supernova remnant (IGelfand fc Gaenslerl 120071 ) . 
Its magnetar nature was confirmed with the detection of radio pulsations with P = 2.096 s and 
P = 2.318xl0 -11 , corr esponding to an infe rred surface dipole mag netic field strength of 2.2 x 10 14 G 



(AXP 1E1547.0 — 5408; ICamilo et al.ll2007l ). An accurate source location was also derived from th e 
radio image, (J2000) R.A.= 15 h 50 m 54 f :ll ± 0?01, decl.= -54°18'23'.'7 ± 0"1 dCamilo et al-lhoOTT ). 



X-ray pulsations at the sam e spin period were later found with a deeper XMM-Newton obser- 
vation (jHalpern et al.ll2008h . No bursts were detected from SGR J1550-5418 until 2008 October 
when both the Swift /Burst Alert Telescope (BAT) and the Fermi /Gamma-ray Burst Monitor 



(GBM ) were triggered by numerous bursts from the source (llsrael et alJl2010i ; Ivon Kienlin et al 
20121 ) . SGR J1550-5418 entered an episode of more active bursting in late 2009 January, and 
no more burst was detected after 2009 April. During these active episodes, several high en- 
ergy instruments, such as the Swift /BAT and X-Ray Telescope (XRT), the Fermi /GBM, the 
Rossi X-ray Timing Explorer (RXTE) /Proportional Counter Array (PCA), and the International 
Gamma-Ray Astrophysics Laboratory (INTEGRAL) /Imager on Board the INT EGRAL Satellite 



(IBIS ) and the SPectromete r on INTEGRAL (SPI) recorded hundreds of bursts (IMereghetti et al 



2009 



Savchenko et al J 12010 ; iKaneko et al.l 120101 : IScholz fc Kaspil 1201 ll ; Ivan der Horst et all 120121 : 



von Kienlin et al.ll2012l ). Following these burst active periods, both the persist ent X-ray emission 



characteristics and the spin-down behavior of the source changed remarkably (lEnoto et al 



Ng et al 



2011 



Bernardini et al. 



2011 



Scholz &; Kaspi 



2011 



Dib et al.ll2012l : iKuiper et al.1 120121 ). 



2010 
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The spectral proper ties of SGR J1550-5418 bursts have been extensively studied using indi 



vidua l instru ments: SPI (IM eregh etti et al.l 120091). IBIS (ISavchenko et al 



20ld ). XRT (jScholz & Kaspill201lh and GBM (|van der Horst et al.l 12012 ). The XRT data cover a 



20ld ). BAT (jlsrael et al 



relat ively narrow en e rgy r ange (0.5 — lOkeV) of the spectrum of a typical SGR burst. Neverthe- 



less, lSj3hj)!z_&Jiaspj (|201ll ) modeled SGR J1550-5418 burst spectra using the XRT data only in the 
energy range of 0.5 — lOkeV with a single power law and found an average photon index of 0.17 
± 0.33. They also reported that there is a slight anti-correlation between the photon index and 
the abso rbed X-ray f l ux. BAT provides a spectral energy coverage for SGR bursts from 15keV to 



150 keV. Ilsrael et al.l (|2010l ) found that the spectra of BAT detected bursts can be well de s cribed by 



a single blackbody function with temperatures ~ lOkeV. Finally, in Ivan der Horst et al.l (|2012l ) we 
derived the spectra for a large set of SGR J1550-5418 bursts detected with GBM in 2009 January 
using several continuum models. We found that in a slightly broader energy range (8 — 200 keV), 
a Comptonization model or the sum of two blackbody functions (BB+BB) can fit equally well the 
SGR J1550-5418 burst spectra. Note that these two models were also used to describe the spectra 



(Feroci et al. 


2004; 


Olive et al. 


2004: 


Israel et al. 


2008: 


Lin et al. 


2011 



In this paper, we combine the spectral data of bursts observed simultaneously with XRT and 
GBM, to investigate their spectral characteristics over a broader energy band (0.5 — 200 keV). In 
particular, we concentrate on the two most plausible representations of SGR burst spectra, namely 
the Comptonization model and the BB+BB. Focusing on bursts with data collected over broader 
spectral bands enhances the chance to discriminate between different spectral models. In Section 
[2j we describe both the XRT and GBM observations of SGR J1550-5418 bursts and the selection 
of their common events sample. We present the data reduction and analysis in Section 12.31 The 
broad band spectral analysis results and their physical interpretation are presented in Sections [3] 
&[U respectively. 



Table 1. XRT observations of SGR J1550-5418 with simultaneous events with GBM. 



Observation ID 


Date 


Start Time 


Exposure a 








(ks) 


00340573000 


2009 January 22 


02:26:22 


6.38 


00340573001 


2009 January 22 


09:18:28 


9.45 


00030956035 


2009 January 30 


17:49:33 


2.97 



a in WT mode 
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Observations and Simultaneous Burst Identification 



2.1. Observations 



T he Swift /XRT is an X -ray imaging spectrometer sensitive to photons in the 0.3—10 keV energy 
range ([Burrows et al.ll2005l ). The telescope is operated either in Photon Counting (PC) mode or in 
Windowed Timing (WT) mode. Both modes provide the same spectral capabilities, however, the 
temporal resolution of the PC mode data is about 2.5 s, too coarse for the very short (< 100ms) 
SGR bursts. Therefore, we only employ here XRT observations performed in WT mode, because 
of its 1.7 ms readout time well suited to the SGR burst durations. XRT monitored the source with 
46 pointed observations in WT mode between 2008 October 01 and 2009 April 30, covering burst 
active episodes of SGR J1550-5418. These XRT observations were densely concentrated around the 
most burst active period (2009 January) with a total exposure time of ~ 175 ks. Three out of these 
46 pointings (listed in Table [TJ included bursts simultaneously detected with GBM, as described 
in detail in the next section. 

The Fermi /GBM monitors the entire sky (excluding the portion occulted by the Earth) in 
the energy range from 8 keV to 1 MeV with twelve Nal detectors and in the 0.2 — 40MeV energy 
band with two BGO detectors. In its trigger mode, GBM records Time- Tagged Event (TTE) data 
with high temporal and spectral re solution of 2 (is and 128 energy channels, respectively. The 
trigger readout lasts for 600 s (see iMeegan et al.l |2009| . for more details of the instrument and 
data types). Using the same burst finding algorithm described in Ivan der Horst et al.l (|2012l ). we 
searched for triggered and un-triggered events from SGR J1550-5418 during the active periods of 
2008 and 2009^ In total we identified 692 bursts out of which, 458 events had TTE data. We 
only used Nal detector TTE data for our spectral analyses, as typical SGR bursts are not detected 
above 200 keV, and last only for a fraction of a second. 



2.2. Identification of Simultaneous Events 

To identify the events observed simultaneously with GBM and XRT, we compared the times 
of the 458 GBM bursts with the time intervals of the XRT observations in WT mode and found 87 
common bursts. Note that some SGR bursts have multi-peaked time profiles and each peak was 
labeled as a burst in our initial un-triggered event search. We used the convention described in 



van der Horst et al.l (|2012l ) to determine whether multiple peaks constituted a single event, namely, 
we requested that the time difference between successive burst peaks was less than a quarter of 
the spin period of SGR J1550-5418 (~ 0.5s). As a result, we obtained 66 SGR J1550-5418 events 
simultaneously observed with GBM and XRT. For each burst, we plotted the burst lightcurves seen 
with GBM and XRT to determine the time interval that includes the main emission episode and is 



x The list of GBM triggered events is available from http://gammaray.nsstc.nasa.gov/gbm/scicnce/magnetars/magnl550triggers.html 
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used as our spectral extraction interval (see Figure [T]). We excluded nine dim bursts that had less 
than forty counts in the XRT data during the burst interval, statistically not enough for spectral 
analysis. Additionally, we e xcluded five very bright bu rsts which saturated the High Speed Science 



Data Bus of GBM (see also Ivan der Horst et al.ll2012l ) 



We checked all common bursts in the XRT data for pile-up. We regenerated the level 2 data by 
including photons in all grades (0 — 15) from level 1 data at the position of SGR J1550-5418, using 
the standard XRT data processing tasks in HE AS OFT. Then we calculated the average readout 
time for each detection ares[§. We found that, besides the five GBM saturated events, in ten other 
bursts the average readout time for central detection areas is smaller than 1.7 ms, the smallest 
readout time for the WT mode. This indicates that these bursts are affected by the photon pile-up 
in the XRT data. Compared to the remaining common events, these bursts have more photons in 
grades higher than 2; as a result an analysis of the data in the good grade range (0 — 2) would 
lack most information and the results would be misleading. Therefore, we also excluded these ten 
piled-up bursts from further investigations. The final outcome of all these filters was a selection of 
42 bursts, observed with both GBM and XRT, that we then used for broadband spectral analysis 
(40 bursts detected on 2009 January 22 and 2 events on 2009 January 30). 



2.3. Data Reduction 

We extracted time integrated spectra for our 42 simultaneously detected bursts. In Figure 
[H we present the burst detected at 02:34:28.194 on 2009 January 22 as an illustrative example of 
spectral integration ranges. We describe below the procedures we followed to generate the GBM 
and XRT spectra. 

For the GBM data, we selected the Nal detectors with a source angle smaller than 60°, and 
without any blockage from the Large Area Telescope (LAT) or other parts of the satellite, such as 
solar panels and radiators. We determined background levels by fitting pre and post burst intervals 
with a first order polynomial using RMFIT v3.4rcl^ and extracted both burst and background 
spectra. We then used grppha to group the extracted source spectrum to include at least 15 source 
counts in each energy bin. We generate the response matrices using GBMRSP vl.9 for each burst. 

For the XRT data, we selected events with grade range of — 2 with xselect, and accumulated 
the source spectra from a 40 pixel long section of the chip centered at the SGR J1550-5418 location 
in the same time intervals as used for the GBM spectra. We extracted background spectra from 
a region of the same size, away from the source. We generated the Ancillary Response Function 
(ARF) file for each burst using xrtmkarf in HEASOFT. In our spectral fitting we used the standard 



2 7 x 1 pixels for XRT WT mode. 

3 R.S. Mallozzi, R.D. Preece, & M.S. Briggs, "RMFIT, A Lightcurve and Spectral Analysis Tool," ©2008 Robert 
D. Preece, University of Alabama in Huntsville, 2008 
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second since 2009 January 22 02:34:28.194 UTC 



1.5 



Fig. 1. — Top panel: Background subtracted lightcurve of a burst detected at 02:34:28.194 on 2009 
January 22 from SGR J1550-5418 with the GBM NaI-6 detector. It is binned with a time resolution 
of 8 ms. Middle panel: XRT lightcurve of the same burst in 8 ms time resolution. Bottom panel: 
hardness ratios (GBM v.s. XRT) of the three sub-intervals indicated with the dotted lines. The 
dot-dashed lines denote the time interval over which the burst spectrum has been accumulated. 
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response file ' swxwt0to2s6-20070901v012.rmf 1 provided in the Swift calibration database. Finally, 
we also grouped the source spectra to include a minimum of 15 counts in each energy bin. We fit 
all spectra using XSPEC vl2.7.0. 



3. Spectral Analysis 



Motivated by our recent results published in lLin et al.l (|201ll ) and Ivan der Horst et al.l (|2012l ). 
we modeled the broadband time-integrated spectra of all 42 common events with the Comptonized 
model (COMPT) and the sum of two blackbody functions (BB+BB). The COMPT model is in a 
single power law shape with a high energy exponential cutoff expressed as: 

/ = Aexp[-E(2 + X)/E pcak ](E/E piv )\ 

where / is the photon flux in photons s _1 cm -2 keV -1 , A is the amplitude with units same as 
/, -Speak is the energy (in keV) at which the spectral distribution function peaks, A is the photon 
index, and E p { v is the pivot energy fixed at 20 keV. The latte r BB+BB model has been commonl y 



used in the context of SGR burst spectra (jFeroci et al.ll2004l : lOlive et al.ll2004l : llsrael et al.ll2008l ). 
In our spectral fits, we fix the multiplicative interstellar absorption term at 3.24 x 10 22 cm -2 , 
since SGR J1550-5418 bursts are short and XRT burst spectra cannot constrain the absorption 
parameter^. In our joint fits, we also include a multiplicative factor to account for the cross- 
calibration between the XRT and GBM spectra. In Figure [2] we show the broadband spectral 
modeling results of the burst of Figure Q] with the COMPT and BB+BB models. We notice that 
a BB+BB model fit to the combined XRT and GBM data results in better residuals than the 
COMPT ones and are also better than the COMPT fits of the GBM data alone (see the two left 
panels in Figure [2]). These results indicate that the broadband joint fits can constrain better the 
physical emission model in SGR bursts as we discuss in Section 14.11 We calculated the average 
reduced x 2 for both models, and obtained 1.03 the BB+BB model and 1.11 for the COMPT model. 
The standard deviation of the reduced x 2 is 0.20 and 0.26 for the BB+BB and the COMPT model, 
respectively. These results further support the BB+BB model against COMPT. 

Table [2] presents the joint fit results with la errors for all 42 bursts. Columns 1 — 3 are the 
burst numbers, start times in UTC, and durations of the time-integrated spectra. The COMPT 
model parameters, i.e., the power law index (A), Ep^, and fit statistics are shown in columns 
4 — 6. Columns 7 — 13 correspond to the temperatures of the two blackbody components, the 
luminosity and size of the emitting area for each blackbody component (assuming a distance to 
SGR J1550-5418 of 5kpc), and the fit statistics of the BB+BB model fits. Since the energy flux 
values obtained with the COMPT model are in agreement with those we get from the BB+BB 



4 We adopt here the value obtained with th e XRT observations of the persistent emission during the burst active 
episode of SGR J1550-5418 bv lScholz fc Kaspil l|201lh 
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50 

Energy (keV) 

Fig. 2. — The spectrum of the SGR J1550-5418 burst shown in Figure [TJ Top two panels show the 
XRT-GBM joint fit spectrum. Bottom two panels show the GBM only spectrum. The left column 
are COMPT model fits and the right column are fits with a BB+BB model. The lower parts in 
each panel show the fit residuals. 
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model, we only present in columns 14 & 15 the energy flux in the GBM energy band (8 — 200 keV) 
and the observed flux in the XRT energy band (0.5 — lOkeV) using the BB+BB fits. 

The COMPT model has one less parameter than the BB+BB function. Since these two models 
are not nested, we cannot employ a A% 2 test to determine which model provides a more adequate 
fit. Therefore, we performed extensive simulations for each of the 42 bursts. We first fit each joint 
spectra with the seed model (i.e., the model with the smallest resulting x 2 value): these were the 
BB+BB model for 33 bursts and the COMPT for the remaining nine bursts). Then, we generated 
10000 simulated spectra based on the seed model of each burst and fitted the simulated spectra 
with both models. For each set of 10000 simulations, we selected fits with well constrained model 
parameters, requiring the errors of E pea k of COMPT and of the two blackbody temperatures to be 
less than 20%. We constructed the distributions of all fit parameters as well as of the fit statistics 
and fitted them with a Gaussian function. The mean values of these distributions for the seed 
model parameters agree with the fit results listed in Table [2j Finally, we calculated the percentage 
of the simulated spectra which result in a smaller x 2 value when fitted with the seed model. This 
percentage, defined as the p— value, reflects the significance of the preference of the seed model at 
a given background and fluctuation level. We conclude that the seed model provides a significantly 
better fit than the other model if p > 0.9. The p— values of all bursts are listed in the last column 
of Table El 

We then grouped the 42 bursts into three categories based on their resulting p— values: the 
BB+BB burst, BB+BB is significantly prefered to COMPT model; the COMPT burst, COMPT 
model is significantly better than BB+BB; and the intermediate group containing bursts for which 
both the BB+BB and COMPT models provide equally acceptable fit results (p < 0.9). We find 
that 31 events are the BB+BB bursts, only one is a COMPT burst (event #18), and 10 are in the 
intermediate group. Figure [3] displays the p— values of all bursts versus their total GBM counts. 
Note that a similar trend is obtained if the p— values were plotted with their corresponding total 
XRT counts. We find that the bright bursts prefer the BB+BB model to the COMPT model. This 
might be an indication of higher opacity, on average, in the more luminous bursts. We explore both 
the statistical and the correlative behavior of all model parameters in the next sections. 

Our joint spectral fits provided the opportunity to investigate the cross-calibration of the XRT 
and GBM instruments with the two spectral models. We determined a multiplicative factor between 
the XRT and GBM detectors for each model. The left panel of Figure U] shows the behavior of this 
factor from the BB+BB model fits as a function of total counts in the brightest GBM detector. The 
same plot for the COMPT model fits is presented in the right panel of Figure HI The instruments 
are perfectly cross-calibrated if the multiplicative factor is equal to 1 (the dotted lines in Figure 
H|). For brighter bursts, the factor from the BB+BB model fit is better constrained, and it does 
not change significantly from burst to burst: its weighted mean value and la error are 1.17 ± 0.05. 
However, the weighted mean of the COMPT model derived factor is 0.58 ± 0.03, much smaller and 
further from 1. We notice that the values obtained from a better fit model are closer to the perfect 
cross-calibration factor of 1. We find that the constant factor is not significantly correlated with 
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Fig. 3. — Plot of the p— value v.s. the total counts in the GBM energy band (8 — 200 keV). The 
black dots are the BB+BB bursts, the open circles are the intermediate group bursts and the five 
point star is the COMPT burst. 
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any spectral parameters. We conclude that the cross-calibration of the XRT and GBM instruments 
works reasonably well within the fluence and energy range of the bursts in our sample. Future 
efforts to better understand the cross-calibration of these two instruments should include bursts 
with wider energy and fluence ranges, and perhaps, different spectral models. 



3.1. Comptonized Model 



We present the distribution of power law indices obtained from the joint fit to the XRT and 
GBM spectra (histograms with thick lines) in the left panel of Figure [SJ We also present in the 
same Figure the distribution of the indices as obtained by fitting the GBM spectra only (histograms 
with thin lines). It is clear that the latter fits yield, on average, lower power law indices (i.e., the 
spectrum is harder), reflecting the overall broadband curvature of the burst spectra. We fit each 
distribution with a normal function and find that the joint broadband fit mean index value is 
-0.58 ± 0.09 (width of 0.43 ± 0.11, dotted lines in Figure U left panel), while the GBM only fit 
mean index is —0.87 ±0.05 (width of 0.42 ±0.06, dashed lines in Figure [5l left panel). We show the 
distribution of the -Epeak values obtained from the joint fit in the right panel of Figure [5j A normal 
function fit to this distribution yields a mean of 45.0 ± 2.1 keV with a width of 10.9 ± 2.2 keV. 
The -Epcak values from joint fit agree with those from GBM data only fit very well, see detailed 
discussion in Section 14.11 

In Figure [H we present a plot of the joint fits -Epeak values v.s. the observed fluence/flux of 
bursts in the 0.5 — 200 keV range. We do not find any an ti-correlation or a broken power law 
trend as see n by fitting the GBM data of SGR J1550-5418 (Ivan der Horst et al.l I2OI2I ) and SGR 
0501±4516 (|Lin et al.l l201ll ). which we attribute to the fact that the 42 common events cover a 
much narrowe r fluence range, about 1/3 of that in the complete GBM burst sample from the same 
active period (jvan der Horst et al.ll2012l ). 



3.2. Two Blackbody Model 

From simulation, we find that the BB±BB model provides a significantly better fit to 31 
common events. We use these BB±BB bursts to investigate the properties of the model parameters. 
In Figure [7] we present the distribution of temperatures of the cool (left panel) and hot (right 
panel) BB components, respectively. We fit a normal function to these distributions and obtain 
a mean value for the cool BB temperatures of 4.4 ± 0.2 keV (width, 0.8 ± 0.1 keV), and for the 
hot BB temperatures of 16.0 ± 0.4 keV (width 2.2 ± 0.4 keV). We also calculated the weighted 
mean values for cool and hot BB temperatures, which are 4.2 ± 0.1 keV and 14.8 ± 0.2 keV, as 
well as their standerd deviations, 0.9 keV and 2.7 keV. These mean temperatures are in agreement 
within uncertainties with t he values obtained by fitting only GBM spectra of SGR J1550-5418 bursts 



(jvan der Horst et al.ll2012l ). and similar to the BB±BB temperature values obtained for the bursts 
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Fig. 4. — Plot of the multiplicative factor values from the BB+BB model (left) and the COMPT 
model (right) fits v.s. the total counts (8 — 200 keV) in the brightest GBM detector. The black 
dots are the BB+BB bursts, the open circles are the intermediate group bursts and the five point 
star is the COMPT burst. The dotted line indicates the factor being equal to 1, and the solid and 
dashed lines display the weighted average of the factor and its la error, respectively. 



12 



10 - 



6 - 



-2.0 



-1.5 



■1.0 -0.5 
hdex A 



0.0 



0.5 



10 



20 



50 



40 



50 
(keV) 



60 



70 



Fig. 5. — Distributions of the COMPT model index (left) and -E pca k (right). The dotted lines show 
the best fit normal functions. The thick histograms in the left panel are the index distribution of 
the XRT-GBM joint fit, while the thinner histograms are the same distribution for the GBM data 
only fit. 
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(Feroci et al. 


2004; 


Olive et al. 


2004; 


Israel et al. 


2008; 


Lin et al. 


2011 



In Figure[8]we demonstrate the anti-correlation between the emission area and the temperature 
of the cool and hot BB components. The Spearman Rank Order Correlation test yields a correlation 
coefficient for the hot BB of —0.79 with a chance probability of 1.40 x 10 -7 . The correlation for the 
cool BB is not as significant as the hot component: its Spearman correlation coefficient is —0.62 
with a probability of 1.75 x 10 -4 . We fit the emission area v.s. temperature for the cool and hot 
BB with power laws, and obtain power law indices of —1.5 ± 1.4 and —4.5 ± 0.9, for the cool and 
hot components, respectively. We also fit the emitting area v.s. the two BB temperatures together 
with a single power law and obtain the best fit power law index of —3.5 ± 0.2 (shown as the solid 
line in Figured]). Note that this value is very close to the theoretical expectation from a single BB 
with fixed luminosity, R 2 cx (ZcT) -4 . The observed departure at the high temperature end from this 
ideal form reflects the higher luminosity present in the hot BB component relative to the cooler 
one; on the average, the hot BB energy is about twice the one emitted from the cool BB. The two 
energies are highly correlated with a Spearman correlation coefficient of 0.88, corresponding to a 
chance probability of 6.45 x 10 -11 ; a power law fit yields an index of 0.99 ± 0.05 (Figured]). 



4. Discussion 

4.1. Impact of the XRT data 

We have used two different instruments, the Swift /XRT and the Fermi /GBM, to perform a 
time-integrated broadband spectral analysis of 42 common events from SGR J1550 -5418. By adding 



the X RT data, we extended the lower energy bound of our earlier spectral analysis (jvan der Horst et al 



20121 ) from ~ lOkeV to ~ 0.5 keV. For most model parameters, our joint fit results agree well with 
the results from fitting the GBM spectra only, as shown in Figure [10] The correlation coefficients 
between all but one (the COMPT power law index) model parameters derived with and without 
the XRT data are larger than 0.94, corresponding to a probability smaller than 3.3 x 10~ 10 . The 
average COMPT index without the XRT data is —0.87 ± .05, consistent within erro rs with the 



mean of ~-0.92 obtained from a much larger burst sample (jvan der Horst et al.l 12012] ). However, 
the inclusion of the XRT data better constrains the COMPT indices, which become harder than 
the ones derived from spectral fits to the GBM data alone, as shown in Figure [5] Therefore, we 
conclude that the COMPT fit to the GBM data only overestimates the emission in the lower energy 
bands. 

Our analysis provides an important diagnostic for the model preference between the COMPT 
and BB+BB models. By adding the XRT data, we find that 31 of 42 bursts are statistically better 
described by BB+BB. This fact, combined with the observed steepening of COMPT model indices 
when excluding XRT data, highlight the generic broad curvature of the SGR J1550-5418 burst 
spectra. We note here that the joint analysis of XRT and GBM spectra is limited to an absorbed 
energy fluence range of 8.3 x 10~ 8 — 1.5 x 10~ 6 erg cm -2 (0.5 — 200 keV), since the brighter events 
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Fig. 6. — Plot of -Epeak v -$- observed energy fluence (left) and average flux (right) between 0.5 — 
200 keV. The black dots are the BB+BB bursts, the open circles are the intermediate group bursts 
and the five point star is the COMPT burst. 
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Fig. 7. — Distributions of the temperatures of the two blackbody components of the BB+BB model. 
The dotted lines show the best fit normal functions. 
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Fig. 8. — Emission area as a function of blackbody temperature for both blackbody components of 
the BB+BB model fit. The hot and cool blackbody components are displayed with triangles and 
circles, respectively. The filled symbols are the BB+BB bursts. The intermediate group bursts are 
shown as open symbols. The solid line indicate the R 2 oc (A;T)~ 3 ' 5 relation, the best fit power law 
function with all emission areas and temperatures for BB+BB group bursts. 
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Fig. 9. — Correlation between the total energy emitted from the hot and cool blackbody compo- 
nents. The filled and open circles present the BB+BB and intermediate group bursts, respectively. 
The best power law fit (E^ot oc -E^ooi) for the BB+BB group is shown as a solid line. The dotted 
line is where the two BB components have equal energy. 



-17- 



would cause pile-up in the XRT data and the dimmer events would not yield high enough statistics 
in the XRT data for a constraining spectral analysis. W e discuss in the nex t session the theoretical 
implications of these results. We also noted the fact that llsrael et al.l (|2008l ) investigated broadband 
spectral properties of a very rare event (the storm) from SGR 1900+14 while our invest igations are 
about much more common typical short bursts. Therefore, our results extend the llsrael et al 
(2008) results to the more common magnetar outbursts. 



We investigated here SGR J1550-5418 burst spectra in a time-integrated manner. As seen in 
the bottom panel of Figure [TJ the hardness ratios of the three parts of the burst, designated by 
the three peaks in the XRT and GBM lightcurves, show a clear hard to soft spectral evolution. 
Detailed time-resolved spectral analysis of SGR bursts would provide important insight to the 
spectral evolution of SGR bursts with time. 



4.2. The BB+BB Model 



The most plausible interpretation of the BB+BB model is the emission originating from two 
hot spots with different temperatures near or on the neutron star surface or in its magnetosphere 
where local thermodynamic equilibria are achieved. It must be emphasized that a BB spectral fit is 
an idealization to that emitted by the physical environment of a real photosphere. Due to possible 
gradients of temperature with optical depth into an evolving region that is approximately in local 
thermodynamic equilibrium, significant distortions from a true blackbody form are predicted in 



SGR photospher ic spectral models (e.g. IUhnerlll994l ; iThompson h Duncan! 1 19951 : iLvubarskv 



Israel et al.l 120081 ) 



2002 



To better understand the BB+BB behavior and uncover its relation with the spin properties 
of SGR J1550-5418 , we investigated the phase characteristics of the 31 BB+BB bursts, as follows: 
we first selected all XRT counts collected during 31 burst intervals and converted their arrival 
times from the Swift mission time to the corresponding time at the Solar system barycenter. We 
then calculated the spin phase for each burst co unt using the appropriate spin ephemeris of epoch 
(MJD) 54854 as reported by bib et al.l (|2012l ) using both RXTE and Swift observations. We 
present the phase distribution of burst counts in the middle panel of Figure [TTJ To ensure that the 
distribution is not dominated by the excessive counts of the brightest bursts, we also calculated the 
probability density for each phase bin, which is the average of the normalized (by total counts) phase 
distributions for all bursts, as shown in the top panel of Figure [TT1 We find that the probability 
distribution of the burst counts is not uniform over the spin phase of SGR J1550-5418 and the 
deviation from the mean probability is significant: we calculate the root-mean-square deviation 
of the phase probability density function from its mean as 0.021 ± 0.001. We also compared the 
phase probability density function to the persistent emissi on phase profile (bottom panel in Figure 
[TT]) obtained using contemporaneous XMM observations ()Dib et al.ll2012l ). The phase probability 
density function is marginally anti-correlated with the persistent emission phase profile in our burst 
sample with the correlation factor of —0.5 corresponding to a chance probability of 3.4 x 10~ 2 . This 
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Fig. 10. — Plots of correlations between model parameters obtained by the joint XRT-GBM fits 
and the parameters obtained by fitting the GBM data only. The black dots are the BB+BB bursts, 
the open circles are the intermediate group bursts and the five point star is the COMPT burst. 
The dotted lines represent the x = y trend. 
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indicates that the burst emission regions on the neutron star surface are not necessarily associated 
with the site persistently emitting in X-rays (typically a BB with a te mperature of 0.5 keV). This 
is in agreement with the crustal fracturing mech anism for SGR bursts ([Thompson &: Duncanlll995l ; 



Braithwaite Spruitll2006l ; iPerna &: Pons!l201ll ) as any portion of the solid crust can fracture if the 



magnetic stress built up is near the threshold to rupture. We also find that the burst probability of 
some spin phases in SGR J1550-5418 is higher. This could be attributed to a non-uniform surface 
magnetic field, with some regions having larger magnetic stresses than others. 



4.3. The COMPT model 



A Comptonization spectrum emerges when low energy photons are repeatedly upscattered by 
the thermal electrons in a corona until the photon energy reaches E ~ kT e . We find that the mean 
value of the Comptonization peak energy, -E pea k is 44.8 keV, which indicates an average temperature 
of the thermal electrons, (T e ) ~ 5.1 x 10 s K. In other words, the average speed of the electrons in 
the corona is ~ 0.4 c, where c is the speed of light. 



The essence of Comptonization spectra is discussed at some length in lLin et al.l (|201ll ). in the 
context of GBM observations of SGR J0501+4516 bursts. While the turnover energy provides a 
diagnostic on the hot electron temperature, the power-law slope below the vF v peak energy defines 
a measure of the op acity in the Comptqnizing region. Specifically, in the simplest theoretical 



Rvbicki &: Lightmanl Il979i ) the index A for a differential photon spectrum 

AkT e /(m e c 2 ) max{r B , r§} 



constructs (e.g., see 

dN/dE oc E x depends only on the magnetic Compton y-parameter ys 
via A 



1/2 — \/9/4 + 4:/ub- Here max{r, r 2 } is the mean number of scatterings per photon by the 
hot electrons, where r is the effective optical depth for scattering, which in our case is modified by 
the strong magnetic field and thus dubbed the magnetic optical depth and denoted by tb • Since we 
require the effective value of tb, one needs to calculate some Rosseland-type mean opacity, averaged 
over photon angles and polarizations and accounting for the effects of the strong magnetic fields 
present, such as anisotro py and polariza tion-mode switching through scattering in non- uniform B 
(see also the discussion in lLin et al.ll201ll ). This index is realized only in the energy range somewhat 
above the soft photon injection energy, Ex, presumed to be surface thermal X-rays, and somewhat 
below the characteristic energy associated with the hot thermal electrons, in this case marked via 
the vF v peak energy .Speak- Moreover, the above relationship for ys requires the scattering to be 
in the Thomson regime, and the mean photon energy to be lower than that of the electrons. 

It is evident that ys S> 1 cases yield the flattest Comptonized spectra with index around 
A ~ — 1. The COMPT fit indices in Table [2] as well as its distribution in the left panel of Figure [5] 
are nearly always harder than this, indicating that repeated Compton upscattering has difficulty in 
generating the observed flat spectra. We note parenth etically, that thi s was also the case for around 
1/3-1/2 the bursts reported for SGR J0501+4516 in lLin et al.1 l|201l] ). The inclusion of XRT data 
in the current work extends the spectral coverage to much lower energies thus enabling us to 
determine the value of A significantly better than in previous works. This results in systematically 
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Fig. 11. — Top: Phase probability density profile of the BB+BB group bursts from 
SGR J1550-5418 ; middle: count phase distribution; bottom: persistent emission pulse profile from 
contemporaneous XMM observations in the 0.5—10 keV band . The horizontal dashed lines in each 
panel represent the mean value of the burst phase probability density (top panel), the burst phase 
counts (middle panel) and the persistent emission count rate (bottom panel). 
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higher values, corresponding to a harder spectral slope, with a mean value of —0.6 ± 0.10 instead 
of —0.92 ± 0.05, which has made this problem worse. The fact that large y# is demanded in this 
fitting protocol would suggest high opacity and strong thermalization might be active in the burst 
emission region. In such cases, the above dependence of ys on T e and tb is not operable, and the 
Compton ization is saturated, and de scribed instead by a modified Wien or modified BB spectrum 



(e.g., see iRybicki Sz Lightmanlll979l ). Moreover, the upscattered power-law photon population is 
generally substantially lower in total number relative to the seed surface X-ray population (i.e., 
it presents itself as an X-ray tail), a broadband spectral shape that is at odds with the spectral 
curvature inferred from the fits here. Hence, there is no strong mandate to prefer a classical, 
unsaturated Comptonization model for the bursts reported here. Accordingly, the spectra might 
naturally be expected to exhibit more truly thermalized character, for example a two-blackbody 



Ulmer 


1994; 


Lvubarskv 


2002) 



imposed by transport within the photosphere (e.g. see lUlmerl 
is what this broad band XRT/GBM analysis has enabled for a magnetar for the first time: the 
clear discrimination between COMPT and BB+BB spectral models. 
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Table 2. Spectral fit result of GBM-XRT common bursts from SGR J1550-5418. 
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56. 


65+ 7 ' 56 
°°-6.05 


56.285/52 


4. 


.731° 


52 
48 


19 


64+1-82 
D4 -1.64 





.4618 


05 

.05 





.881° 


08 
.08 


7.151 3 ; 4 ? 


26 


04: 


33: 


:29. 


.134 





.168 


01+ ' 42 


43 


I 7+2.87 

II -2.78 


42.645/31 


1 


7118 


72 
51 


11. 


64l°; u l 


1 


1218 


(>(> 
57 


1 


.9818 


14 
.14 


1025.33+gg 39 


27 


04: 


33: 


:35. 


.750 





.480 


—0 66+°' 23 
U.DO_ 22 


42 


n +2.95 
1 -2.77 


68.651/59 


4. 


1318 


71 

69 


14. 


.24ll; 4 8 


0. 


.5818 


08 

"■> 


1 


•1818 


10 
11 


15 05+15.02 ' 


28 


04: 


33: 


50 


.294 





.560 


_n 74+0 28 
u - ' -0.26 


42. 


5 o+4.82 
08 -4.11 


44.346/51 


5 


.3818 


59 
5fi 


19 


oq+3.04 
za -2.58 





.5618 


07 
07 





•6718 


09 

.09 


c 91 +2.32 
°- zi -1.52 


29 


04: 


34: 


12 


.430 





128 


fl oq+0.35 
u - z ' :5 -0.32 


47. 


oq + 3.38 

uc> -3.10 


36.254/32 


4. 


•24tS 


15 
93 


13. 







77+° 


23 
18 


2 


■5518 


24 
29 


I8.5ll? 7 ;f 2 


30 


04: 


34: 


14. 


502 





392 


-0 24+°' 47 
u - z *-0.42 


56 


04+7.41 
' 84 -5.86 


28.03/22 


3 


.2018 


99 
89 


14. 







1918 


07 

.06 





■9218 


10 
.10 


13.98l 46 4 46 


31 


04: 


34: 


20. 


.622 


1 


.184 


-0 36+ 10 


42 


40+1. 06 
41 -1.03 


219.531/158 


4. 


.6218 


24 
23 


14. 


fi6 +0.53 
DD -0.50 


1 


.2318 


08 

07 


2 


.3318 


09 
09 


21.131 3 ; 94 


32 


04: 


39: 


23. 


636 





.432 


-0 64+°' 25 
U,D -0.23 


59. 


-7+6.60 
°'-5.29 


41.121/50 


6 


.3311 


52 
fi3 


19 


4 g+4.29 
^°-3.22 





4918 


17 
15 





.9618 


14 
.16 


9 07+2. 08 



21+ '" 6 
u - zi -0.05 

12+ ' 14 
u - lz -0.07 

0.46+ ' 14 
0.02 
0.03 



0.35 
0.09 



0.11 

+0.04 
0.02 
+0.04 
0.02 
+0.17 
-0.12 
0.04 
0.03 
+0.22 
0.09 
+0.15 
0.12 
+0.17 
-0.14 
0.20 
0.15 



0.16 
0.49 
0.64 
0.54 

1 1 1+0.26 
1-11-0.23 
n 1O+0.05 
u.io_ 04 

0.77+ ' 24 



0.31 



0.22 
+0.12 

0.09 



0.25 
0.08 
0.26 
0.35 
0.10 



c+0.19 



-0.12 
0.06 
0.04 
+0.12 
0.09 
+0.12 
0.10 
+0.05 
-0.04 

09+ ' 04 
u - uy -0.03 

in+°. 2 ° 

u - lu -0.07 

16+ 003 
U-lD_o.03 
5+0.23 
-0.17 
-+0.02 
-0.01 
-+0.22 
*-0.19 
,+0.11 
-0.08 

04+°'° 3 
U - U4 _0.02 
+0.24 
0.20 
+0.06 
-0.05 
0.07 
0.06 
+0.07 
0.03 



0.48Z 
0.05^ 
0.84j 
0.22 4 



0.53 
0.15 
0.39 
0.05 



158.980/131 

29.642/38 

69.707/67 

46.837/49 

26.868/27 

41.331/31 

22.212/28 

12.402/18 

55.320/50 

26.889/33 

80.685/75 

114.304/100 

29.603/26 

43.926/43 

66.447/61 

26.934/23 

35.231/27 

53.494/45 

23.320/26 

29.538/45 

42.803/36 

37.277/27 

134.705/130 

62.295/59 

47.499/51 

34.018/30 

67.225/58 

39.108/50 

34.438/31 

25.48/21 

165.813/157 

49.488/49 



7.51 
1.28 
3.06 
1.85 
1.11 
1.04 
1.06 
0.71 
2.14 
1.40 
3.27 
5.24 
1.04 
1.57 
2.65 
0.96 
1.03 
1.69 
1.21 
1.96 
1.60 
1.01 
7.73 
2.49 
2.37 
1.18 
2.55 
2.13 
1.32 
3.50 
12.9 
2.01 



+0 
-0 


25 
.28 


0. 


.9218 


04 
07 





.9996 


+0 
'-0 


12 
18 





.2218 


01 

09 





.9766 
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.10 
12 
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.1718 
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.09 
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28 
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.4218 
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.3518 
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Table 2 — Continued 



COMPT COMPT COMPT cool BB hot BB cool BB hot BB cool BB hot BB BB+BB 

^peak 



# T start a dt A E peak x 2 /dof kT kT L 39 b L 39 b R 2 R 2 x 2 /dof F GBM C F^ nT d 



(UTC) (s) (keV) (keV) (kcV) (km 2 ) (km 2 ) 



33 04: 

34 04: 

35 05: 

36 05: 

37 06: 

38 07: 

39 10: 

40 10: 

41 19: 

42 19: 



39:29.908 
39:36.628 
49:17.289 
49:46.105 
06:41.596 
26:28.960 
49:29.944 
52:11.888 
29:41.189 
29:42.941 



0.184 
0.336 
0.536 
0.200 
1.000 
1.672 
0.456 
0.248 
0.432 
0.120 



+0.19 
0.18 
-I -.0+0.23 
-1.18_ 21 

n on+0.20 
u " 5z -0.20 
+0.49 
0.44 
+0.19 
0.19 
+ 0.10 
0.10 
+0.22 
0.21 
+0.48 
0.46 

-0 54+ - 12 
U.O4_ 12 

+0.44 

0.39 



-0.71 



-1.08] 
0.171 
-0.56^ 
-1.251 
-0.5ll 
0.541 

-oW 



45.33 
24.37 



41.63±i;ig 
36.93lg \% 

+3!05 
2.82 
+4.45 
6.21 

25.49tMf 
59.821^20 

54.941^1 
66.48t°;!» 



77.969/67 

50.047/53 

91.248/76 

41.745/27 

82.334/69 

317.180/180 

53.961/50 

9.511/16 

204.166/160 

27.940/33 



4.07 
4.16 
4.47 
1.91 
3.84 
4.23 
3.90 
3.99 
4.19 
5.42 



+0.35 
0.33 
+0.37 
-0.34 
+0.68 
0.61 
+0.57 
0.43 
+0.57 
0.59 
+ 0.18 
0.17 
+0.39 
0.36 
+ 1.28 
1.53 
+0.37 
0.34 
+ 1.76 
1.51 



14.05t8£ 

15 21+ 1 ' 61 
+ 1.00 
-0.86 



9. or 

12.00+?,;^ 

1 ft 1 q+0-58 
lO.iy_0.55 

20.38tl'gf 
14.26l«« 
16.85±g;|| 



1.42 
0.87 
0.57 



+0.14 
0.14 
0.07 
0.08 
+0.09 
0.09 



1.13 



19.29 



+ 2.1 



+ 1.01 
0.41 

75+o.iQ 

1.07 
0.78 
0.38 
1.07 
0.45 



0.12 
+0.04 
0.04 
0.07 
0.07 
+0.12 
0.13 
+0.07 
0.07 
+0.16 
0.13 



2.87 
1.10 
1.15 
1.10 
0.63 
1.70 
1.28 
0.57 
2 
1.75 



0.18 
0.19 
+0.10 
0.10 
+0.10 
0.11 
+0.11 
0.12 
+0.12 
0.11 
+0.06 
0.06 
+0.11 
0.11 
+0.14 
-0.15 
0.10 
0.10 
+0.18 
0.20 



40.24t£06 
22.80±|;|g 

ll.ost™ 

666.351 2 ™ 

26.46™ 

+4.42 
3.72 

26.27l^ 3 6 f 

+61.30 
-7.26 

26.87±3, 1 g2 3 

+9.94 
2.49 



0.57 



+0.15 



25.99] 



11. t 



4.05Z 



-0.12 
,+0.06 
' -0.04 
,+0.09 
-0.06 
l.ou_ 4g 

U - Z4 -0.15 

n iq+o-03 
u.iy_ .03 

.+0.02 
-0.02 
+0.18 
-0.08 
-0.04 
-0.04 
+0.05 
0.04 



0.13Z 
0.171 



0.061 
O.lll 
0.28 
0.10 



64.717/66 

35.201/52 

90.915/75 

38.411/26 

70.145/68 

180.111/179 

37.806/49 

10.539/15 

189.646/159 

24.662/32 



QO + U.IJM 

z - J °-o.io 

1 95+° 10 
1 - ao -0.12 
9 oq+0.15 

+0.05 
0.12 
+0.18 
0.22 
+ 0.28 
0.32 
+0.16 
0.16 
+0.09 
0.11 
+0.13 
0.13 
+0.05 
0.06 



0.87Z 

3.791 

13.851 

2.761 

0.691 

5.231 

0.85^ 



0.31 
0.35 
0.49 
0.11 
0.87 
1.88 
0.31 
0.20 
0.32 
0.12 



+0.04 
0.04 
+0.04 
-0.05 
+0.04 
0.05 
+0.00 
-0.09 
+0.04 
0.07 
+0.08 
0.08 
+0.03 
0.04 
+0.01 
0.07 
+0.03 
0.03 
+0.01 
0.04 



0.9585 
0.9943 
0.9868 
0.9913 
0.9990 
1.0000 
0.9998 
0.2544 
1.0000 
0.9313 



a Events 01 — 40 were happened on 2009 January 22, while bursts 41 and 42 were on 2009 January 30. 

b The luminocity of Blackbody components in units of 10 39 ergs -1 

c Butst energy fluence in GBM band (8 — 200 keV) in units of 10 — 7 ergcm~ 2 

d Butst energy fluence in XRT band (0.5 — lOkeV) in units of 10 -7 ergcm~ 2 



to 



